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SIGNIFICANCE  AND  EXPLANATION 


We  consider  the  motion  of  a  viscous  fluid  filling  the  whole  space  R*, 
governed  by  the  classical  Navier-Stokes  equations  (1).  Existence  of  global 
(in  time)  regular  solutions  for  that  system  of  non-linear  partial  differential 
equations,  is  still  an  open  problem.  From  either  the  mathematical  and  the 
physical  point  of  view,  an  interesting  property  is  the  stability  (or  not)  of 
the  (eventual)  global  regular  solutions.  Here,  we  assume  that  v.j(t,x)  is  a 
solution,  with  initial  data  a^x).  For  small  perturbations  of  a^,  we  want 
the  solution  v^(t,x)  being  slightly  perturbed,  too.  Due  to  viscosity,  it  is 
even  expected  that  the  perturbed  solution  v2(t,x)  approaches  the  unperturbed 
one,  as  time  goes  to  +«.  This  is  just  the  result  proved  in  this  paper.  To 
measure  the  distance  between  v^(t,x)  and  v2(t,x),  at  each  time  t, 
suitable  norms  are  introduced  (L^-norms). 

For  fluids  filling  a  bounded  vessel,  exponential  decay  of  the  above 


distance,  is  expected.  Such  a  strong  result  is  not  reasonable,  for  fluids 
filling  the  entire  space./  In  this  case,  we  prove  that  the  I^-distance 
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Other  notations  are  standard,  or  will  be  introduced  in  the  sequel.  Let 
p  >  3  and  a^  e  L1  n  LP+2,  div  aj  =  0.  Assume  that  there  exists  a  global 

oo  p+2 

solution  e  L  (0,+®;L^  )  of  problem  (1),  with  initial  velocity  a^  and 

pressure  tt^.  This  solution  is  strong  and  unique.1  We  prove  the  following 
stability  result: 


i  p 

Theorem  A  Assume  that  the  above  conditions  hold  and  let  a2  e  L  n  L  ,  be 
such  that  div  a2  =  0.  Then,  there  exists  a  positive  constant  y0  such  that 
if 

(2)  la i  -  a2 1 p  <  Yq  # 

there  exists  a  unique  solution  v2  £  C( [0,+“>) ;lP )  of  (1)  with  initial  data 
a2*  This  solution  verifies  the  estimate 

(3)  | v  1  ( t )  -  v2(t)|  <  C(1  +  t)"3/4  . 

1  2 

The  constants  yQ  and  C  depend  on  p,  on  the  L  and  L  norms  of  the 

CO  p+2 , 

initial  data  a1  and  a2,  and  on  the  L  (0,+<»;Lr  )  norm  of  v.,.  In 

particular,  by  considering  initial  data  a2  such  that 

a2  1 1  <  lai  I  i  +  k1  '  I a2  ! 2  <  l31  I2  +  k2  ' 

where  k1  and  k2  are  any  positive  constants,  y0  and  C  depend  only  on 

k.,  k2,  and  on  the  norms  !aili»  | a  1 1 2  and  flv,B  +2  of  the 

!  L  (0,+»;LP  ) 


unperturbed  solution  v^. 


1 

We  refer  the  reader  to  the  results  proved  in  [2).  See  also  [5],  and 
references  there. 
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The  local  existence  and  uniqueness  of  a  strongly  continuous  solution 
v2(t)  with  values  in  L2  n  iP  is  well  known.  The  bound  (3)  guarantees  the 
global  existence  of  v2(t). 

The  proof  of  theorem  A  follows  the  method  introduced  in  reference  [1]  in 
order  to  study  the  asymptotic  behaviour  of  the  solutions  of  system  (1). 


Proof  of  theorem  A.  The  difference  w  *  v2  -  v^  satisfies  the  following 
system: 


w'  +  (v2*V)w  +  (w^VJv^  -  Aw  +  VP  =  0 


div  w  =  0 


w(t=0  "  a 


in  R  , 


where  P  =  tt2  -  ir1  and  a  =  -  a^.  Multiply  both  sides  of  equation  (4)  by 

|w|P-2w  and  integrate  over  R^.  After  suitable  integrations  by  parts  we 


obtain 


-  |w|P  +  N  (w)  +  4  /  |  7  |w|p/2|2  = 

p  dt  p  p  2  1  1  11  ' 

P 


-  /  (v2*V)w*|w|P  2w  -  /( w»V)v^  •  | w | 2w  -  /VP»|w|P  2w 


The  first  term  on  the  right-hand  side  is  zero  since  v2  is  divergence  free. 
By  integrating  by  parts  the  other  two  terms  we  get 


-  af  lwlp  +  <  (p  -  D  J  | w | P“ 1 | Vw | | v 1 | 


(p  -  2)  /  |P| |w|p_2|Vwj  . 

Consider  the  first  integral  on  the  right-hand  side  of  (5).  By  Holder's  and 
Young's  inequalities  one  has 
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(p  -  1)  /  MP-'IVwMvJ  <  (p  -  1)N  (w)  (/MPjv^2)1/2 


g  N  (w)  +  2 (p  -  I)2  /  |w|P  [vj2 


<  1  N  (w)  +  2(p  -  1)2|w|P  2  |Vl 


On  the  other  hand,  one  can  prove  the  following  inequality  (see  [1],  equation 


(1.14)): 


3  p(p-1) 

|w|p+2  <  c  Vw)  P+2  |w!P  ?+2  * 


hence,  from  (6),  one  obtains 


(p  -  1)  /  |v|P'1|Vw|lv1|  <  l  Np(w)  +  c|w|P|Vl|p+2 


p-1 


Consider  now  the  second  integral  on  the  right-hand  side  of  (5).  By  Holder's 
and  Young's  inequalities  we  have 

(p  -  2)  /  |p|  |w|P'2|<7wj  <  (p  -  2)(/|p|2|w|P“2)  1/2  Np(w)V2 


<  2(p  -  2 ) 2  /  | P| 2 |w| P-2  +  |  Np(w) 


<  2(p  -  2)2|P|^2l«l^2  ♦  g  »p<“l  • 
2 


To  estimate  P,  we  take  the  divergence  of  (4)  and  obtain 


-  *P  =  J.  3jkwl(V1]  +  V23)  =  J.  WIT  ^i(2vi3  + 

ifD  i  D  id  i  3 


From  the  Calder<5n-Zy ground  inequality  one  has 


|p|Pi2  * c  T.i'i<2vi +  > 

2  2 


i*TV 
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then,  by  Holder's  inequality. 


1*^2  <C<“lp+2<|v,|^2*  W^2 


)  . 


By  introducing  this  last  inequality  in  (9)  and  by  using  inequality  (7),  we  get 
(p  -  2)  /|p||w|P  2 1 Vw|  <  cMj*2(lVl!p+2  +  lwlp+2)  +  I  NP(W) 


(10) 


3  p(p-1)  3 

<  C  N  (v)P+2|wf  P+2  KJ2  +  C  N  (w)P|w|P"1+  l  N  (w) 


2  JB+2.  P^P-1) 

1  _  p-1  P-3 

<  -  Np(w)  +  c|  w |P  |v  ^  +  C  |w|p 


Hence,  by  using  inequalities  (8)  and  (10),  we  get  from  (5) 


(ID 


2(p+2)  p(p-1) 

-  ■—  | w |P  +  r  N  (w)  <  c|w|P |v  |  -  +  c|w | 

p  dt  1  'p  2  p  P  1 'p+2  1  'p 


On  the  other  hand,  one  can  prove  the  following  Sobolev  type  inequality 
(see[1],  equation  (3.2)) 

Np(w)  >  C|w|Pp  , 
and,  by  interpolation,  one  has 


4  3(p-2 ) 

Mp  <  M?'2  lwi3pP"2 


Hence 


(12) 


N  (w)  >  c|w|”8|w|*+8  , 


'P 


where  8  »  4p/3(p-2).  On  the  other  hand,  by  a  result  proved  in  reference  [31 

2 

(see  also  [4])  on  the  L  -decay  of  the  solutions  of  the  Navier-Stokes 


equations,  one  has 


(13)  |w(t) | 2  <  | v 1 ( t ) | 2  +  | v2( t) | 2  <  C  (t  +  1)~3/4  , 

1  2 

where  C  depends  only  on  the  L  and  L  -norms  of  the  initial  velocities. 
Then,  from  (12)  and  (13),  we  obtain 

(14)  N  (w)  >  C  (t  +  1)38/4|wIP+f?  . 

P  P 


Hence,  from  (11)  and  (14),  we  have 


1  d£  |w|p  +  c(t  +  i)3S/4iwip+B  <  c|w!P!vi|P+r1  +  c>wip p_3  * 


2(p+2 ) 


>(p-1 ) 


p 1  1 'p+2 


from  which  we  obtain,  since  v^  is  bounded  in  L P+2,  uniformly  in  time. 


w  +  C. 


t(t  +  1)36/4|w|8  +  1  <  C 2 1 w |  +  C3|w 


1+B+Y 

P 


where  Y  =  2p  /3(p  -  2)(p  -  3).  In  (15),  C1  depends  on  p  and  on  the  L 

2 

and  L  -norms  of  the  initial  velocities,  C2  depends  on  p  and 

v  1  ^  +2  '  C3  dePen<Js  only  on  p.  Consider  now  the  corresponding 

L  (0,+°»;LP  ) 


o  ,  d. e  . 


y’(t)  +  C 1  ( t  +  1)3B/4(y(t))R  +  1  =  C2y(t)  +  y(t))1+8n 


y (  0 )  =  |  a| 


We  prove  now  that,  if  | a | ^  is  sufficiently  small,  then 


y(t)  <  C(t  +  1) 


By  comparison  theorems  for  o.d.e.  it  will  follow  that 


| w( t ) |  <  y(t);  hence  (3)  holds.  Let  tg  e  ]0,+°°[  be  such  that 


c  +  c 

*„  >  (-V-*) 4/36  -  1 


WVWOVWr 


’  Ml 


**  v  y 


By  the  continuous  dependence  on  the  initial  data  of  the  solution  of  (16),  one 
can  find  yQ  >  0  (depending  on  p  and  on  C^,  i  =  1,2,3)  sufficiently  small 
such  that,  if  lalp  <  Y Qr  then  y(t)  <  1  for  each  t  e  [0,tQ].  Moreover,  if 
there  exists  t  >  t^  such  that  y(t)  =  1,  then  from  (16)  and  (17)  one  has 

y* (t)  =  -  c^t  +  i)36/4  +  c2  +  c3 


<  -  C^tQ  +  1)36/4  +  c2  +  c3  <  0  . 


This  implies  y(t)  <  1,  for  any  t  >  0.  From  (16),  we  then  obtain 
y(t)  <  z(t)  for  any  t  >  0,  where  z(t)  is  the  solution  of 


z'(t)  +  C^t  +  1)38/4(z(t))B+1  =  (C2  +  C3)z(t)  , 


(18) 


z(0)  =  |a|p  . 


Equation  (18)  is  of  Bernoulli  type  and  its  solution  is  given  by 


(C  -k:  )t  t  8(C?+C^)s  , 

z(t)  -  e  [  | cx |  +  f  e  BC^(s  +  1)  ds] 

P  0 


1 

e 


Hence 


*  t,38/4  < 


—oo  f 

y°  *  Bc’  J 


t  0(C2  +  C3)s  36/4 

e  ( s  +  1 )  ds 


By  the  i'Hopital  theorem  one  easily  shows  that  the  right-hand-side  of  the 
above  inequality  converges  to  (C2  +  C3)/C^  as  t  -*•  +“•  Since  it  is  equal  to 
yQ  for  t  =  0,  it  follows  that  it  is  bounded,  in  the  interval  [0,+®)  by  a 
constant  C  (which  depends  only  on  C^,  C2,  C3  and  p). 
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Ja^  -  a2jp  <  y0,  it  corresponds  a  (unique)  global 
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